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ABSTRACT 

This document presents some results obtained with a digital com- 
puter program written at Goddard Space Flight Center to obtain electro- 
magnetic fields scattered by perfectly reflecting surfaces. For purposes 
of illustration a paraboloidal reflector is illuminated at radio frequencies 
in the simulation for both receiving and transmitting modes of operation. 
Fields are computed in the Fresnel and Fraunhofer regions. A dual- 
reflector system (Cassegrain) is also simulated for the transmitting case, 
and fields are computed in the Fraunhofer region. 

Appended results include derivations which show that the vector 
Kirchhoff-Kottler formulation has an equivalent form requiring only inci- 
dent magnetic fields as a driving function. Satisfaction of the radiation 
conditions at infinity by the equivalent form is demonstrated by a conver- 
sion from Cartesian to spherical vector operators. A subsequent develop- 
ment presents the formulation by which Fresnel or Fraunhofer patterns 
are obtainable for dual-reflector systems. A discussion of the time- 
average Poynting vector is also appended. 
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GLOSSARY OF NOTATION 


Symbol 
E(x', y', z 


')H(x',y', z') 
fj., e, a 


j. i 

1* H l’ E 2’ H 2 

'P 

V«/> 
d s 

d -C 

n 

e., h. 

P 

v 

0 

k 

\ 


e 

4> 

<t> 

D 


<P> 


R 

F 


Meaning 

backscattered electric and magnetic fields at observer 
position (x 1 , y',z') 

constitutive parameters: magnetic permeability, inductive 
capacity, electric conductivity 

angular frequency 

imaginary operator = sj -1 

electric and magnetic illumination distribution fields 

solution to the wave equation 

vector gradient of \fj 

differential area 

differential arc -length 

unit normal to a surface 

electric and magnetic polarization vectors 

vector from a feed point to a surface 

a vector used to generate e . and h. for a spherical wave 

order, origin 

wave number = 2 n/k 

wavelength 

radial distance in spherical system 
polar angle in spherical system 
azimuthal angle in spherical system 
unit vectors in spherical system 
aperture diameter 
time-average Poynting vector 
distance from origin O 
focal length of a paraboloid 

vii 



A feed displacement in focal-plane 
K beam deviation factor 
9 beam squint angle 
K electric sheet current 
, electric surface and line charges 

a, i radial and angular variables for surfaces in parametric 
form 
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THE CALCULATION OF ELECTROMAGNETIC FIELDS 
IN THE FRESNEL AND FRAUNHOFER REGIONS 
USING NUMERICAL INTEGRATION METHODS 


INTRODUCTION 

The purpose of this document is (1) to present some computed results ob- 
tained with the Kirchhoff-Kottler formulation in the Fresnel and Fraunhofer 
regions, (2) to discuss equivalent scattering formulations, (3) to introduce a 
simplification in the procedure for obtaining Fresnel and Fraunhofer patterns 
for dual-reflector systems, and (4) to derive the expression used for the time- 
average Poynting vector in the digital computer program. 

The results obtained here are based on the contributions of many persons. 

It is impossible to acknowledge all of them properly. Most of the fundamental 
ideas were well-developed many years ago, but the advent of the fast digital 
computer has, within the past 5 years, made possible the numerical evaluation 
of existing formulations for practical geometries, feeds etc. Since this docu- 
ment contains representative results in Fresnel as well as Fraunhofer regions, 
some of the language and concepts of optics are introduced. The bright diffrac- 
tion disc and rings of Airy, for example, are especially significant in the recep- 
tion mode of operation (image space). A brief historical outline is included. 
Most of the dates and comments in this introduction were taken from textbooks, 
journal articles, etc. 


History of Diffraction 


1452-1519 

da Vinci: 

first reference to diffraction 

1601-1665 

Fermat: 

principle of least time which implied rectilinear 
propagation, independent rays, law of reflection, 
law of refraction 

1621 

Snell: 

sometimes credited with laws of reflection and 
refraction 

1663 

Gregory: 

"Gregorian" dual-reflector system 

1663 

Cassegrain: 

"Cassegrainian" dual -reflector system 

1665 

Grimaldi: 

first accurate descriptions of diffraction 

1690 

Huygens : 

first proponent of wave theory without periodicity 
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1642-1727 

Newton: 

corpuscular theory (received acceptance for 100 years) 
to the exclusion of Huygens wave theory 

1718 

Maraldi: 

found the "Arago" spot in the umbral region 

1817 

Young: 

suggested transverse rather than longitudinal wave 
motion and the interference principle 

1818-1827 

Fresnel: 

extended Huygens' theory by adding periodicity in 
time and space to the wavefronts, but several argu- 
ments had no physical significance. Accounted for 
rectilinear propagation as \->0, reflection, refraction, 
and diffraction effects. Predicted bright spot in disc 
umbra 

1818 

Poisson: 

reviewed Fresnel's prize memoir, and refuted conclu- 
sion of bright spot in shadow of disc 

1818 

Arago: 

verified Fresnel's bright spot by experiment 

1827 

Fraunhofer: 

diffraction formulations, measurements 

1824-1844 

Hamilton: 

geometrical optics (ignored Huygens' and Newton's 
formulation; indifferent to physical interpretation 
and interference concepts) 

1834 

Airy: 

calculation of diffraction disc and rings 

1849 

Stokes: 

limited version of general treatment later made by 
Kirchhoff 

1865 

Maxwell: 

wrote the fundamental equations of electromagnetic 
theory V x E = - BB/Bt and V x H = J + BD/3 t 

1882 

Kirchhoff: 

wrote the Kirchhoff-Helmholtz Integral Theorem, 


1883 

1888 


Poynting : 
Maggi: 


0 (x', y ' , z') 


1 

4 Tf 



B \p 

B n 


d s. 


Kirchhoff started with the scaler wave equation in 
three dimensions, and deduced a formulation which 
in many respects embodies the basic ideas of the 
Huygens -Fresnel principle. 


theorem on energy flow 

provided additional insight on Kirchhoff theory, in- 
cluding a transformation for an unclosed surface 
with rim F. 
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1893 

Heaviside: 

studied the energy flow problem and physical inter- 
pretation of formulations 

1909 

Debye : 

analytical work on scattering 

1909 

Reiche: 

investigated axial non- linear phase shift 

1911 

Sommerfeld: 

rigorous analytical solutions, used boundary value 
approach 

1917 

Rubinowicz: 

preceded Kottler in studies of unclosed surfaces with 
rim 

1923 

Kottler: 

showed that Kirehhoff's solution to the black-screen 
diffraction problem was a first approximation to a 
boundary value problem (actually a solution to a 
"saltus" problem). Kottler annexed a contour inte- 
gral, for unclosed surfaces, to the Kirchoff formulation. 

1939 

Stratton and 
Chu: 

further analysis of Kottler's contour integral and its 
relationship to Maxwell's equations 

1957 

Keller: 

presented an extension of the theory of geometrical 
optics 

1967 

Sancer: 

showed that the vector Kirchhoff equations can be 
derived using the free-space dyadic Green's function 
and the appropriate Green's theorem. The Kottler 
contour integral was obtained without appeal to physi- 
cal intuition, and was inherently contained in the 
derivation rather than appended at the close (the 
historical development. 


The preceding partial list of contributors is indicative of the broad field of 
analysis and research associated with the subject of diffraction. In this document 
the Kirchhoff-Kottler formulation has been stressed with a view toward obtaining 
meaningful solutions for practical feed-reflector configurations. Exploration of 
the bounds of validity of the selected theory is predicated on objective compari- 
sons with other analyses and direct physical measurements where practicable. 
The interested reader may wish to review the mathematical and physical rea- 
soning which underlie the Kirchhoff-Kottler theory of diffraction. Problems 
associated with the definition of a black screen, the introduction of polarization 
in a physical 3-space, and finite surfaces are discussed in detail in numerous 
textbooks 1 and will not be considered here. The extrapolation from the original 

1 Ref. 1, Chapter ||. 

Ref. 2, pp. 460-470 
Ref. 3, Chapter 5 
Ref. 4, Chapter VIII 
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black-screen derivations to present day applications, such as those involving 
multimode waveguide feeds in dual-reflector systems, provides a sharp contrast. 
A review of the analytic difficulties encountered in the theory, the departure 
from the assumption of continuity required by the underlying Green's theorem, 
and the additional compromises due to numerical evaluation should be weighed 
against nearly a century of useful predictions resulting from the Kirchhoff theory. 
Barakat 1 , in paraphrasing a remark of Poincar£, said "The theoreticians believe 
in the Kirchhoff theory because they hold it to be an experimental fact, while the 
experimentalists think it to be a mathematical theorem." The statement is a 
tribute to the theory considering that scatterer size, observer range, and radii 
of curvature of both scatterer and wavefront must be large in terms of a wave- 
length. Further, the effect of mutual coupling on the scatterer, and between 
source and scatterer, are traditionally neglected altogether. 


FORMULATION OF THE SCATTERED FIELDS 

The formulation for the scattered electric and magnetic fields from per- 
fectly conducting surfaces is taken as 2 


E (x', y', z') = ^_L_ _L 

j w e 4 77 


[(n xH.) - V]Vi//ds-jaj/i JU (n x H. 

4 77 J 1 


) 0 d s 


H (x', y', z') = - -J_ (n x H.) x V 0 d s, 

4 77 4 


where 


p-ikr / i 

0 = ^— , V0 = - (j k + i)01 r 


Previously, the equivalent electric -field formulation 


E (x', y', z') = 


- — — ~ V 0 H. • d % - _L [j o)/z(n x H.) 0 + (n- E.) V0] d s 

J ui M 77 1 477) 1 1 

J C • / S 


1 Ref. 5 Chapter 9, pp. 6-16 

Appendix A provides background for this formulation from several sources. 
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was employed. 1 The latter formulation for E(x',y',z') could be termed the 
"historical" approach since the Kottler contour integral appears explicitly as 
an amendment to the Kirchhoff formulation. It has been shown that this line- 
charge integral can be obtained along with the Kirchhoff formulation without 
appeal to physical intuition via a dyadic Green's function, 2 and further, the two 
forms for E(x',y',z') yield identical results. 


From a programmer's viewpoint the organization of the problem is compli- 
cated by the "historical" formulation for the electric field E(x',y',z'). It is noted 
that a mixture of contour and surface integrals requires separate quantizing or 
sampling intervals for numerical evaluation. In addition, the surface-charge 
integral appearing in the "historical" formulation involves the incident electric 
field, so that the polarization vector e. must also be developed and retained in 
the numerical evaluation. This should be contrasted with the formulation used 
herein, which proceeds simply with surface integrals, the incident magnetic field, 
and the associated magnetic polarization vector h. for the incident wave. 


e. 

I 


P J X (V X p t ) 

I V (vxPj)l' 


h. 

1 


Pj X V 

IPj X vj 


Unless there is a specific reason for distinguishing the Kottler line -charge 
integral and Kirchhoff surface-charge integral contributions in a problem, their 
combined influence can be computed via the surface integral 


1 1 

j co e 4 7T 


[(n x Hj) • V! V \p d s 
•4i 


for the total bound electric charge distribution. It is easy to verify that the 
formulation for E(x',y',z') and H(x',y',z') satisfy the radiation condition at in- 
finity. A slight complication is encountered with the operator del, written 


V = 


3 

3 x 


— + k — 
3 y 3 z 


in Cartesian coordinates. 


1 Ref. 2, page 469 

2 Ref. 6 
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The radial component, at infinity, of the integral 


- )OJ fX 


4 Ti 


f (n X H x ) 

s i 


can be set down immediately as 


The integral 


- j 


1 

co !± 

4 co jx 



L( n x H x ) 0 - l r ] l r cl s 


_1 1_ 

j oo e 4 77 



[( n x H x ) • Vj V \jj d s 


is entirely radial at infinity, as is anticipated for charges, however, manipulation 
of the integrand in terms of the Cartesian operator V is undesirable since 

is radial . 1 An interpretation of (n x H x ) • V in spherical coordinates is, therefore, 
needed to establish the vanishing of the radial part of the electric field at infinity. 

A vector dot product in orthogonal curvilinear coordinates can be defined 


P • q = P u Q u + P v q v + P w q w 

and the dot product (n x H x ) -V should follow directly except for the fact that 
p -> (n x H x ) is a vector, and q “* V is an operator. Furthermore, the form of V 
in spherical coordinates depends on the operation to be performed, which is 
unlike the Cartesian case where V is the same for grad, div, curl, cross-product, 
and dot-product. 


1 The numerical evaluation of [(n x H x ) 'V] ■ V './/also needs to be considered. A conversion of 
V i/rfrom spherical to Cartesian coordinates was used previously in the program, posing no 
problems, but in the present formulation 


(hx H,) -V =K 
V x 


3 

3 x 


+ K 

y 



_ 3 _ 
3 z 


is a differential operator , and V0 is a numeri ca I result for each differential area ds. 
2 Ref. 7, page 155. 
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The form of V to be used is suggested by the quantity on which scalar form 
[ (n x ) • V] operates. Here 

j h + I) 0 l r 

is a function of (r) alone, and is a vector . The spherical basis vector l r is 
invariant with radial changes. From orthogonal curvilinear coordinate relation- 
ships, 1 the gradient of a scalar (£) in spherical coordinate is 



Vf = l Hi 

a 3r r 0 


U i 

d0 9 


1 

7 ^ 


i. 


and 


a = 


d r 


( r K) 


= i 


Then 



r 



r s i n 8. 


[(n x Hj)- V] V0 = 


(nx Hj) • l r 


jT 

d* r_ 


v </» = [(n 


s,)-i r ] 

3 r 2 


= l r {- k 2 t2ib + -p [(nx H t ) • i r ] 

which shows (l/r) in orders O 1 O 2 , O 3 since 0 = e"i kr /r. The highest order 
of the term (l/r) in the formulation for Efx'.y'.z') containing the contour integral 
is O 2 


1 Ref. 7, page 152 
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Vanishing the radial fields at infinity is now made very simple. 1 Summing 
the radial parts of 


E (x', y', z') = — — [(n x H.) • V] V 0 d s - j u/i J— (nxH)^ds 

K ) e 4 7T J 4 77 J s 

b l 1 

at infinity, terms of O 2 , O 3 in (l/r) are neglected. Remaining terms of O 1 in 
(l/r) annihilate. 


1 

4 TT 


k ^ - j 0 )fj. (n X Hj)- l f 


j a) e 


If = o 


and 


c = f A. =. 


vl/x e 


an identity, is recovered. 

The magnetic field is inherently simpler due to the absence of bound mag- 
netic charges. 


H (x', 



r 

(n x Hj) x V \}j d s. 

4 , 


This field is transverse since V0 is radial at infinity and the cross-product of 
two vectors is orthogonal to the vectors entering into that cross-product. 


1 The numerical evaluation can also be carried out now since the integrand [(n x H-|) *V J V <p 

has been rewritten as l f (-k^ + 2 j k/r + 2/r^)0 which is convertibile to Cartesian components 

and free of troublesome operators. V \p-- (j k + l/r) >/> 1 is assumed to be reduced to the 
Cartesian or “working notation” of the program already. 
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NOT REPRODUCIBLE 


Figure 1. NASA Tracking Station, Rosman, N.C. 
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SINGLE-REFLECTOR SYSTEMS 


The vector Kirchhoff-Kottler diffraction theory is applied here to the study 
of paraboloidal systems, initially, then to the conics, and subsequently to dis- 
torted conics and more general surfaces exhibiting aberrations of various kinds. 
Results presented in this document are principally for the paraboloid in both 
transmit and receive modes of operation. 1 The paraboloid is singular in its 
focussing properties, has found wide application in optics and microwave tech- 
nology, and extensive data exists for it in the literature. Particularly interesting 
are those paraboloids with small f numbers (f = F/D < 0.5) as these are encountered 
often in radio -frequency applications and are somewhat more difficult to analyze 2 
than the large F/D structures. 

A typical RAD HAZ environment exists in the installation at the NASA track- 
ing station (Rosman, N. C.) shown in Fig. 1. Some years ago (circa 1964) a study 
of the possible radiation hazard existing on one of the access roads was under- 
taken. The only methods available at that time were awkard and empirical. 3 In 
general, the field-strength and energy-density levels at a set of points in the 
vicinity of a high-gain paraboloid energized by a high-power transmitter were 
sought. Without trying to re-create the precise set of conditions at Rosman, an 
attempt will be made to illustrate how the Diffraction Program developed at 
Goddard Space Flight Center might be used, not only in a RAD-HAZ application, 
but also in a more general way to design feeds for scanned beams, explore other 
reflector geometries, and capture energy from the Airy -type diffractions rings, 
etc. 


To illustrate the RAD-HAZ application: Fig. 2 and Fig. 3 show the trans- 
mitted electric field versus angle 9 as the radius of the observer is reduced 
from 100 miles through 2D 2 A, 1D 2 A, and D 2 /2>v . Coordinates in the field 
explored by these four arcs are representative of the access road mentioned 
previously. The object is to determine amplitude, and to a lesser degree polar- 
ization and phase in a RAD-HAZ study. It can be seen that the field pattern shape 
is only slightly degraded in the Fraunhofer region between 100 miles and 2D 2 /\ . 


1 The results presented in this document formed the basis of a presentation give at the 1971 
USNC-URSI Spring Meeting, Washington, D. C. under Commission I, Near-Field Measurements 
and Radiation Hazards Probing (RAD HAZ). 

2 Ref. 1 1, page 485 
Ref. 12, page 353 

Ref. 13, see section on antennas (power density multiplication) 
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F = 12.99 FT, D =30.0 FT, f = 1.11520 GHZ 




Figure 3. Electric Field vs. Angle 6 
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Some null-filling is in evidence, the rr -radian phase jumps at pattern nulls have 
become more gradual transitions at pattern minima and signal level is nearly 
proportional to l/R. Within the 2D 2 /A arc the pattern is unique for each radius, 
null-filling is pronounced, beamwidth is increased, and the signal level no longer 
varies as the inverse first-power of distance. 

As the paraboloid is approached more closely, the moving spherical triad of 
basis vectors conforms poorly to the transmitted wave which, in the focal plane 
(R = 0), resembles a plane wave more than a spherical wave. For this reason 
the Cartesian basis is used to resolve the fields. In addition, the time-averaged 
Poynting vector is now displayed. 1 See Fig. 4. The region considered is centered 
about the axial point D 2 /2\ from the focus of the paraboloid. To the far-left of 
the area shown by Fig. 4, the more intense time-average Poynting vectors < P> 
form a nearly parallel bundle whose diameter is roughly that of the paraboloid. 

To the far-right of the area shown, the vectors < P> diverge and appear to come 
from origin O, which is also the focal -point F of the paraboloid. That is, the 
set of vectors <P > in the Fraunhofer region appear to flow radially from a point 
origin, but their intensity is modulated in accordance with the envelope shown by 
Fig. 2. Fig. 4 is interesting because it shows that transition region where the 
Poynting vectors < P > have not quite aligned themselves with the lines tracing 
back to origin O. No attempt was made to scale the length of the power flow 
vectors < P> , but it was noted that those vectors which departed considerably 
from the lines intersecting at O were also reduced in intensity. 

A measured result, obtained for a region between the focal point and a point 
approximately 17 X further away from the transmitting paraboloid, is shown in 
the text by Collin and Zucker. 2 It is noted that although the discussion has been 
limited to paraboloids with focal-point excitation, the program is written in a 
general way to accept other surfaces, including distorted conics, etc. and there 
is no restriction on the number of feeds, method of excitation, position of feeds 
etc., other than those restrictions inherently imposed by the Kirchhoff theory 
itself. 

At this point it is convenient to invoke reciprocity and study some of the 
fields obtained from reflectors in the receiving mode of operation. The received 
fields are less intense than the transmitted fields, above, and do not constitute a 
RAD-HAZ problem in the preceding context, however, an element of RAD-HAZis 
still present since, for the case of paraboloids, spheres and other reflectors, 
energy is being concentrated in a relatively small region. From a standpoint of 
feed design, the received fields provide a means of matching feeds to fields (fre- 
quently in a conjugate sense). The literature of optics abounds in both experi- 
mental and theoretical studies of the image space . 3 As before, the paraboloid 

1 Appendix B 

2 Ref. 14, page 50 

2 Ref. 4, pp. 435-480 
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is utilized here for convenience and because its properties are well-known. An 
objective of the program, however, is to obtain insight for those scattering sys- 
tems that are difficult to treat analytically and for which only meager data can 
be found in the literature. 

The constant field intensity plots (isophotes) shown in Fig. 5 are for the 
focal-region of an ideal paraboloid for the axial reception of a plane-wave. The 
source in this case is + o the right, the scatterer to the left, and the converging 
spherical wave is progressing from left to right. This convention will be re- 
tained in all of the subsequent examples. Results are symmetrical about the 
y-z plane, but asymmetrical about the focal plane. The F/D of this system is 
only 0.433, therefore, energy flows into the image space at a steep angle. Fig. 6 
shows the Airy disc and rings in the x-y on focal plane. 

The phase of the electric field is shown in Fig. 7 for the preceding set of 
isophotes (Fig. 5) over a restricted portion of the z-axis (±1A. re O). Power- 
flow in this region is orthogonal to these wave fronts. It is also apparent that 
there is no rational definition of phase for axial nulls, and zero-intensity rings 
of the Airy pattern. Energy flow is around such singular points. 

A plot of the direction of the time-average Poynting vector <P> for the 
focal region being discussed is given as Fig. 8. The crossed lines are the bounds 
of the geometrical shadow. Airy discs, bright rings, and dark rings can be seen 
without difficulty. No attempt was made here to scale the length of vector < P > 
in accordance with its magnitude, however, this capability has now been auto- 
mated in a Calcomp plot program. The turbulent regions to the left and right of 
the geometrical focus correspond to the limits of the depth of field at which the 
intensity of field vanishes. It is noted that the time-average Poynting vectors 
do not cross the system axis, which is distinct from a ray-optics presentation, 
but rather diverge as if they originated from F for large z. 

Especially interesting is the region surrounding a null. The null between 
the Airy bright disc and the first bright ring is shown in Fig. 9. A circulation 
or vortex-like flow can be seen about the null and an examination of the magni- 
tude of <P> shows that the intensity of the time-average Poynting vector dimin- 
ishes in a continuous manner as the null in the vortex is approached. The 
existence of this counter-flow of energy has been reported previously, and was 
computed by Minnett and Thomas 1 who used axial -wave theory to investigate the 
characteristics of fields in the focal region of a paraboloidal reflector. 


1 Ref. 15 
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FOCAL PLANE 


Figure 8. Direction of the time-average Pointing vector (axial reception) 



f=1 .11520 GHz, 0=30.0 ft., F=12.99 ft. 


INS 

o 


X X. X X 

XXX X 

X X -x X 
\ \ \ X 

\ \ \ \ 

* \ \ \ 

\ \ \ \ 

\ \ \ \ 

\ \ \ \ 

\ \ \ \ 
\ \ \ \ 

-0.2X 


x 

\ 

* 

I 

I 

\ 

X 


X 


\ 

\ 

\ 

\ 

\ 


X 


\ 



\ 


x, x ~x ~- 

X / / / 
X X X ^ _/ t f ft 

w u /_A \ u 

I I / i y ^ w\\ 

l t f '' v\\U 

w\\^~ 

\ X x X 
X X X "■ 



X 


jr 

X 

y 

y 

X 


y 

y 

X 

y 

y 

y 

y 

y 

y 

y 

y 

/ 

y 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

t 

t 

t 

/ 

/ 

t 

/ 

t 

f 

t 

t 

t 

t 

f 

\ 

t 

t 

t 

t 

t 

t 

t 

t 

f 

t 

t 

t 

t 

t 

t 

t 

f 

t 

t 

f 

/ 

/ 

/ 

f 

/ 

t 

t 

y 

y 


y 

/ 

/ 

/ 


- 0.8X 


- 0.7X 


- 0.6X 


-0.1X Z=0X 0.1X 


0.2X 


FOCAL PLANE 


Figure 9. 


Direction of the time-average Pointing vector in the vicinity of a null (axial reception) 


Additional information concerning the calculation and measurement of am- 
plitude and phase in the image space (of a microwave lens) has been given by 
Famell, , Bachynski, and Bekefi, 2 who applied scalar theory. Emphasis is placed 
on the fact that vector Kirchhoff-Kottler theory was used to obtain results pre- 
sented in this document and that a numerical approach has been stressed to 
avoid many of the mathematical difficulties encountered in analysis, particularly 
for those cases where the diffracting system was complicated by intricate feeds, 
or illuminating sources, and non-conic or distorted surfaces. The latter are 
accommodated by the program without additional effort on the part of the pro- 
gram user, who simply modifies the input parameters, and the penalty in in- 
creased computer running-time (CPU time) is usually moderate. 

A very simple illustration of the flexibility of the program is afforded by 
displacing the illuminating source so that a plane-wave (effectively) arrives with 
an arbitrary angle of inclination relative to the system axis (z-axis). The simu- 
lation produces this change for the source by calling Euler-angles in a rigid-body 
transformation and the components of a displacement vector. An inclination 
angle of 4-degrees was chosen, together with the same input parameters used 
previously for axial reception, based on the anticipated displacement of the in- 
tense field region associated with a focussing paraboloid. The equation 


s in 9 

S 


~ KA 


F 


relates beam squint-angle (? s through be am -deviation factor K 3 , feed displacement 
A, and focal length F. Under reciprocity then, for F/D = 0.433, 


s i n 4° = ( -89) (A) 

(12.99) 

and the energy can be expected to concentrate near 

A* 1.0 2 ft. => 1. 14 \ (S) 1.11520 GHz, 
on the opposite side of the z-axis. 


1 Ref. 16 
Ref. 17 

2 Ref. 18 

3 Ref. 3, page 488 
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Figure 13. Direction of the time-average vector (paraxial reception) 


A displaced and distorted version of the isophote field for axial reception 
is shown in Fig. 10. It appears that the peak intensity occurs just to the left of 
focal plane (toward the reflector). Solid lines at 4° are shown through F and F' 
as a reference for the displaced high intensity focal region. Fig. 11 shows the 
displaced Airy disc and rings in the x-y or focal plane. The wavefronts for the 
Fig. 10 computations are shown in Fig. 12. As before, the time-average Poynting 
vector was computed for the paraxial case. See Fig. 13. The Airy bright disc is 
displaced, together with the bright and dark rings. As with previous examples, 
power flow in the vicinity of dark rings and limits of the depth of field, which 
tends to be orthogonal to the general trend from left to right, is accompanied by 
large reductions in magnitude of the power density. The vectors <P> were not 
scaled at this writing because of the large amount of data reduction involved for 
this asymmetric example. 


DUAL-REFLECTOR SYSTEMS 

An initial effort to obtain the fields scattered from dual-reflector systems 
was undertaken by computing the Fraunhofer fields from a Cassegrain configu- 
ration. See Appendix C. Several idealizations were made, including omission of 
multiple reflections between y and y 2 , mutual coupling of the sheet currents, 
and even first-order blockage effects. The latter effects were taken into account 
previously in single -reflector systems by deleting the integration over the umbral 
region. A polygonal approximation to the shadow boundary was specified for 
each observer position. 

The radiation pattern scattered by the main (parabolic) reflector was ob- 
tained as follows. A prime source 3 t = S, eos N 0 was selected with N = 35.0 so 
that the feed taper plus space taper of the illumination at the edge of the (hyper- 
bolic) subreflector equaled -9.0 db. The electric sheet current on the hyper- 
boloid 7 j was then computed with contiguous increments of area not larger than 
0.3 k by 0.3X. Each value of sheet current K on the paraboloid y was then 
obtained by integrating over all of K i on y i , 2 A sampling grid not 2 larger than 
3*0 k by 3.0 A. was used on the main reflector. Finally the distribution K was 
integrated in the Fraunhofer region for each observer position (r\ 9‘ , The 
distributions Kj and K 2 are invariant with observer position and were there- 
fore stored in the computer. For the reflector and feed geometries employed 
here, a four-fold symmetry exists for and K ^ however, this feature was 
not exploited even though cpu time can be conserved in such situations. 
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The results for the dual-reflector system, computed by applying Kirchhoff 
theory redundantly to the surfaces, are shown as Fig. 14. The degraded second- 
null is probably duetto the sampling intervals chosen for integration of the dis- 
tributions Kj and K 2 . The 8 domain for Fig. 14 was adequately sampled 
(A 8 =0.01 deg.) to detect the null, if it developed in the computations. This can 
be seen by an inspection of the phase characteristic which is also degraded 
in the vicinity of the second null. It is a more or less general rule that sampling 
intervals for sheet current distribution must be smaller for large values of 8 . 

The topics of sheet-current-distribution sampling and cpu time were also 
explored from the standpoint of stationary phase. For sufficiently high fre- 
quencies the geometric optics solution should be approached and, therefore, 
one might search for stationarity with the idea of reducing cpu time. The radial 
and azimuthal parameters, o- and £ , used to generate the hyperboloid can be 
suitably restricted on y 1 for a point on y 2 at which K is to be evaluated. By 
this method only a portion of y l is integrated - hopefully all of the stationary 
part, to the exclusion of the self-annihilating or non-stationary terms. The 
search for stationary phase regions was restricted by control over the radial 
variable (a) alone in the initial attempt, but will be repeated with control over 
£ , the angular variable, subsequently, to form increments of area. A typical 
increment of area on the hyperbolic subreflector will then be integrated by 
specifying a . , a , £ . , £ , and the observer point on -y (for that region 

of Kj on y ) will be determined by geometrical optics. Preliminary results for 
one particular set of surface and feed parameters indicated that stationary con- 
tributions ranged over most of y 2 when the observer position was in the vicinity 
of the vertex of y 2 . When the observer position was near the edge or limbs of 
y 2 > however, only an annular region of y x appeared stationary. 

A separate study of Cassegrain and Gregorian subsystems consisting of 
prime feed and subreflector (y ) was carried out to explore the geo- 
metrical optics conclusion that such a subsystem is effectively a virtual point 
source in the limit as A. - 0. The feed 3 is at the conjugate focus F* of y x 
and the virtual source should appear at F of y x , which also coincides with F 
of y 2 in the fully assembled system. Since the parameters used in the studies 
were not those of Fig. 14, graphical results are not presented here. For a 
hyperboloid diameter of 24 k and an ellipsoid diameter of 32 X it was found that 
the radius or range of observation relative to focus F could be reduced until, 
for wide angles, the distribution K x or y x was no longer optically visible. 

That is, the Cassegrain or Gregorian subsystems behaved as virtual point 
sources situated at F. Amplitude variations on a spherical locus of observation 
departed from a monotonic descent by about ±1 db typically. Phase variations 
on the same spherical locus of observation varied by approximately ±3 electrical 
degrees from the phase of a point source situated at F. In conclusion, as the 
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radius of observation is reduced (R -* 0) and the frequency is increased (k - 0) 
the observer does not lie in the near-zone (2D 2 /A.) for the conics studied here 
since the scattering approximates that from a virtual source and D is effectively 
zero. 

Preliminary calculations have been made for the "receive" mode of opera- 
tion for dual reflectors, applying the theory of Appendix C with near-field inte- 
grals present. The Airy disc and rings have been observed in the conjugate focal 
plane for a coarse sampling (LI). A mapping of the focal region is planned for the 
example shown in Fig. 14, taking into account the domain of interest in the x-y 
and x-z planes for a magnification factor M = (e + l)/(e - 1). The optimum 
sampling (LI LI 2 ) of the distributions and K 2 is not known at this time, 
and several computer runs will be required to establish a reasonable tradeoff 
between cpu time and accuracy of results. The magnified width of the Airy disc 
is a factor here since LI must usually be decreased as the observer moves off 
of the system axis. Over-sampling of Kj and K 2 eventually establishes the 
"stability" of the solution with respect to LI, observer position, etc. 


SUMMARY 

The electric and magnetic fields, their associated phases, and the time- 
average Poynting vector were computed in the Fresnel and Fraunhofer regions 
for several well-known scatterers. Transmit and receive modes were discussed 
although, from a program standpoint, there is no logical distinction. A dual- 
reflector system was considered, and it appears that even multiple reflectors 
might be treated by the methods presented, however, the blockage (obscuration) 
problem requires additional effort. 

Present plans call for extensive application of the existing program to ground- 
based and spacecraft antennas. Measured prime-feed data will be injected into 
the simulation and distorted antennas of various types will be treated. A variety 
of methods for obtaining directive gain and spillover will be compared and docu- 
mented. Long-term plans include studies of multibeam polarization-diversity 
and frequency-diversity systems, shaped reflector and Schwarzschild configura- 
tions, parabolic cylinder antennas, spherical reflectors, and doubly curved 
surfaces. This document is an interim report. 
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APPENDIX A 


THE FORMULATION OF THE SCATTERED FIELDS 
E(x',y',z') H(x',y',z') BY EQUIVALENT METHODS 


The following formulation of Stratton 1 was used, initially, to obtain fields, 
time-average Poynting vectors, isophotes and wavefronts 


E (x‘, y', z') = - _ — L_ _L_ (h 

j m e 4 7? J 1 4 


bj 

Js. 


[j (Ofj. (n X Hj) I p + (n-E,) V0] ds 


H (x\ y', z') 


nfj (n 

J Sj 


x Hj)x Vi/<d s 


^ = > v ^ - (j k + 1 ) ^ i r . 


An immediate conversion to the formulation used subsequently can be 
effected by making use of a relationship in one of Stratton's proofs. 2 

J [( n x H) - V) V <// d s = - j> V^Hj • d l - j w e J (n-E)Vi//ds 


From the above, the equivalent form of the elctric field is 


(nx H 1 )i//d s. 

The operator Vy^ of Stratton can be avoided by the use of dyadic notation. 3 It is 
noted that the equivalent form for the scattered electric field was obtained by 
Stratton in the course of demonstrating that the "historical" form of E(x',y',z’) 
satisfies Maxwell's equations. 


E (x', y', z') = 


1 1_ 

j co e 4 77 


[( n x H) V] Vi//d s - j w/i 


4 TT 


^ Ref. 2, page 469 

2 Ref. 2, page 470 

3 Ref. 8, pages 27-28 


33 


The equivalent form of the elctric field is also recognized in the general 
dyadic derivation of Sancer. 1 2 Attention's called to the fact that Sancer con- 
cludes with the forms of E(x',y',z') and H(x\y\z') of Stratton, and does not dis- 
cuss the more compact_formulation which disposes of the contour integral, the 
incident electric field E £ , and its polarization vector e. . 

In a recent text by Rusch and Potter, 2 the equivalent compact form for 
E(x',y',z') is presented together with a discussion of open surfaces. Sign vari- 
ations appear between the formulations of Stratton, Sancer, Rusch and Potter, 
and Silver. 3 These appear to be due to conventions pertaining to the Green's 
function i/> and its gradient, V \p. In practice, for perfectly reflecting surfaces, 
the E(x',y',z') form has two integrals which must annihilate the radial fields and 
the H(x',y',z') form has only one integral. Signs, therefore, are not a vital issue 
here even though easy comparisons are occasionally frustrated by the appearance 
of what seems to be a spurious minus sign. 

The text by Silver (1949) presents the "historical" form of E(x',y',z') 4 and 
the compact equivalent form 5 as a specialization from volume to surface inte- 
gration. 6 The "historical" form evolves when a single surface integral is equated 
to the sum of a surface-change integral plus a line-change integral. One might 
speculate on the consequences of a complete transformation of E(x',y',z') to line- 
integral form 7 for open and closed surfaces. 


1 Ref. 6, page 142, equation (2.17) 

2 Ref. 9, page 46-47 

3 Ref. 3, page 160 

4 Ref. 3, page 160 

5 Ref. 3, page 132 

6 Ref. 3, page 87 

^ Ref. 1, page 79 
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APPENDIX B 


THE TIME -AVERAGE POYNTING VECTOR 
The time-average Poynting vector is defined 1,2 as 


<P)> = Re E (x', y', z') x Re H (a', y', z') = I. Re E (x', y' 


z') x H* (x', y', z') 


It is noted that either the real part or the imaginary part of a complex solution 
(such as E or H) may be chosen at the conclusion of a calculation to represent 
a physical state. For squares and products, however, the real part of the result 
is not equal to the product of the real parts, of E" and H, say. Details of the 
derivation of < P >can be found in several textbooks . 3 

The "working notation" of the digital computer program is the Cartesian 
vector basis, and the complex aspects of the solution are retained by separating 
real and imaginary components. After integration of the reflector surface the 
computer stores 


E xR’ 

E xr 

E yR’ 

£ yr 

E z r- 

E r 

Z I 

H xR> 

H xr 

H y*’ 

H y i’ 

H z r- 

H y 

z I 


the real and imaginary vector components of the electric and magnetic fields. 
Then the time-average Poynting vector 


<P> = I Re E (x\ y', z') x H* (x' , y' , z ' ) 


-Re [i (E y H ; - E z Hj) + j (E z H* - E x H* z ) + k (E x H] 


E y h;)j 


= i Re [ i {(E yR + j E yI ) (H zR - j H zI ) - (E zR + j E zl ) (H yR - j HyI )} 


1 Ref. 2, pp. |131-137 

2 Ref. 10, pp. 370-377 

3 Ref. 2, p. 136 
Ref. 3, p. 70 
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?r> 


+ 1 •> + k {. ...}] 

= — Re 
2 

1 ^ E yR H zR + E yi H,I - E zR H yR - E zI H yI ) + j (E yI H zR - E yR H zI 

“ E zl ^yR + E zR ^yl^ + j U ’ ' ’ ) + J ( )] + k [(...) + j (....)] 

~ 2 ^ ^ E yR + E y i H 2I - E zR H yR - E zl Hyj) 

+ J ( £ zr Hxr + E zi H xl - E XR H zR - E XI H ZI ) 

( E xR ^yR + E xl HyI " E yR H xR - E yI H xl )] 

gives the time-average power flow in watts/meter 2 . 
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APPENDIX C 

FRAUNHOFER FIELDS FOR DUAL-REFLECTOR SYSTEMS 1 


OBSERVER (X', Y z') 



Kj, q qp 
1 1 


(*i» y j- *j) 

Figure C-l. Abstract Dual-reflector system 

Consider a dual-reflector system as in Fig. C-l. The feed is designated 
3j , main- reflector y 2 , and subreflector r y 1 , Associated sheet currents, sur- 
face changes, line charges, and normals as given by K, q s , q^, and n appropri- 
ately subscripted. The fields due to 3 X are proportional to E , H , the fields 
due to are proportional to E 2 , H 2 , and the scattered fields due to y 2 are 
designated E(x',y T ,z'),.H(x',y',z'). 

1 Due to Raymons Miezis of Computer Applications, Inc. 
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Assume initially that every point of y 2 lies deep in the Fresnel region of y 1 , 
and that the observer is in the Fresnel region of y 7 . The general formulation is 
examined to determine what integrals must be evaluated. 


E (x', y', z') = — 


j CO £ 4 77 


[( n 2 x H 2 ) • V] V0 - j 


COfl 


J S, 


4 77 


( n 2 x H 2 ) 4> d s 


H (x', y', z') = - J_ f (n x H, 

4 77 * 

*'S. 


) x V d s. 


Now only H 2 is required above, and 


h 2 ( x 2> y 2 - z 2 ) = - 


4 77 


(n x x Hj) x V 0 d s. 


Since the observer is assumed to be in the Fraunhofer region of y 2 here, 
the transverse fields are obtained from 


E (“, 6\ $•) = - j M fi _L 

4 IT 



(n 2 x H 2 ) \p d s 


projected onto the basis vectors l g , and 1 .. H(co,<9', 4>‘) is most economically 
obtained from the auxiliary relationship 


[ E (oo, d\ <j >‘) | -_ 

|h («>, d', 0')| z ° 

In this formulation only H 2 is required and only two integrals, (x 2 , y 2 , z 2 ) and 
E (co, Q', cf. >') are evaluated since the observer is in the Fraunhofer region, com- 
pared with four integrals for the general case, observer in the Fresnel zone. 

The advantages of the compact equivalent form of the electric field can now 
be seen. If E(x’,y',z') were written in the "historical" form 


E(x', y', z') = — _ 


1 


J coe 


1 

477 


Vl /'H 2 - d^ - _L [(j cvfi(n 2 x H 2 ) 0 + n 2 - E 2 ) V i/r] ds 

J C ^ vc 
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for the Fresnel region of y 2 , E 2 as well as H 2 would be required for the evalu- 
ation. But 


) V \p\ d s , 

which is complicated by the appearance of a contour integral and involves a mix- 
ture of electric and magnetic fields in the integrands. 


0 2 - y 2’ ^2^ 


j coe 4rr 


Vi/zHj • dl l -_L 


[j ^(r^x 


39 


